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Abstract 

A mechanical analogy is used to analyze the interaction between the magnetic field, electric 
current and deformation of interfaces in liquid metal batteries. It is found that, during charging 
or discharging, a sufficiently large battery is prone to instabilities of two types. One is similar to 
the metal pad instability known for aluminum reduction cells. Another type is new. It is related 
to the destabilizing effect of the Lorentz force formed by the azimuthal magnetic field induced by 
the base current and the current perturbations caused by the local variations of the thickness of 
the electrolyte layer. 

PACS numbers: 47.20.-k,47.20.Ma, 47.65.-d 
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FIG. 1: A scheme of a liquid metal battery. Three liquid layers B, E, and A fill a cavity, which may 
be of cylindrical, 3D rectangular, or other shape. During the charging and discharging processes 
the uniform vertical electrical current of density Jo = Joe z driven by electrons and positive ions of 
metal A is imposed. It generates a purely or approximately azimuthal magnetic field (not shown). 
Other components of the magnetic field can be generated by electrical currents in neighboring 
batteries and supply lines. 

I. INTRODUCTION 

The work presented in this paper is motivated by the efforts to develop the liquid metal 
battery, a device for short-term stationary energy storage. Small laboratory prototypes have 
already been shown to work and demonstrated potential for higher efficiency and longer 
operational life than the traditional solid-electrode batteries (see, e.g., fl,Q). The key 
question now appears to be whether larger, more efficient, and commercially viable devices 
based on the same principle can be designed. 

A simplified scheme of the battery is shown in Fig. [D It is a vessel filled with three liquid 
layers: liquid anode A made of a light metal (e.g., Na, Li, or Mg) at the top, layer E of a 
molten salt electrolyte in the middle, and liquid cathode B containing a mixture of a heavy 
metal (e.g., Bi, Sb, or PbSb) and the compound between the heavy and light metals at the 
bottom. The electrolyte is chosen so that it is immiscible with the liquids on either side and 
conductive to the positive ions of the light metal. The system is stably stratified, with the 
density of the light metal being about two times smaller than the density of the electrolyte 
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and many times smaller than the density of the heavy metal. The top and bottom walls of 
the vessel serve as current collectors, while the sidewalls are electrically insulating. 

The electric energy stored in the battery is the difference in the Gibbs free energy between 
the light metal in its free state and in compound with the heavy metal. The processes of 
charging or discharging correspond to, respectively, the electrochemical reduction of the light 
metal from the compound and forming the compound. The reactions occur in liquid state, 
at the interfaces between the electrolyte and the metal, and in the presence of strong (about 
1 A/cm 2 ) electrical currents flowing in the vertical direction. 

An inspection of the scheme in Fig. |T] suggests that the operation of a large-scale battery 
will differ significantly from that of a small laboratory prototype. The reason is the hydro- 
dynamic instabilities which will appear and become stronger at larger size. The result of the 
instabilities will be fluid flows in all the three layers with the potential implications for the 
battery operation that can be both positive (enhanced mixing of reactants) and negative 
(spatial and temporal non-uniformity of reaction rates and, in the worst case, deformation 
of the interfaces so strong that it leads to rupture of the electrolyte layer and disruption of 
the operation). This issue has been a subject of close attention recently. Several mecha¬ 
nisms of the instability have been identified, such as the Tayler instability (sen, e.g., Bfl), 
electrovortex instability 4], and thermal convection caused by bottom heating [6] or internal 
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Joule heating of the electrolyte [71. It has been confirmed that the instabilities are active 
in batteries of even modest size (of radius about 20 cm in the case of Tayler instability and 
as small as a few cm in the case of internal heating convection). Further investigations are 
needed to fully understand the instabilities and reveal their effect on the battery’s operation. 

Yet another likely instability mechanism, which has not been considered before, is ad¬ 
dressed in this paper. It has magneto-electro-hydrodynamic nature and is related to the 
fact that, during the charging or discharging processes, strong current passes through liquid 
layers of vastly different electric conductivities. The conductivity oe of the electrolyte is 
about four orders of magnitude lower than the conductivities a a and as of both metals. This 
means that even a small deformation of the electrolyte-metal interface, i.e., a small variation 
the local thickness of the electrolyte causes a strong variation of the local resistance and, 
thus, significant changes in the distribution of the electric currents within the battery. In 
this paper, we explore the possibility that the Lorentz forces resulting from the interaction 
of the electric current perturbations and the magnetic field act on the liquids in such a way 
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that the deformation of the interface is enhanced. 

On the level of basic physics, the concept of such an instability is not new. A similar 
mechanism has been found in the Hall-Heronlt aluminum reduction cells, where it is called 
the ‘metal pad instability’. A reduction cell is a horizontally large (about 3 by 10 m) 
and shallow (about 20-40 cm) rectangular bath filled with molten aluminum at the bottom 
and molten salt electrolyte with aluminum oxide dissolved in it at the top. The ratio of 
electric conductivities between the metal and the electrolyte is about the same as in liquid 
metal batteries. Electric current of density about 0.1 A/cm 2 flows predominantly vertically 
through the two layers causing the desired effect of electrochemical reduction of aluminum 
from its oxide. 

For many decades, the aluminum industry faced the major problem of the metal pad 
instability that developed in the form of growing sloshing waves at the aluminum-electrolyte 
interface when the distance between the top electrode and the interface was too small or 
the current density was too high. If allowed to evolve, the instability led to short circuit 
between the aluminum and the top electrode, in which case the operation of the cell had 
to be stopped. Keeping the thickness of the electrolyte layer above the threshold resolved 
the problem, but at the cost of substantial energy losses to the waste Jonle heating of the 
electrolyte. Remarkably, proneness to the instability varied among the cells of the same 
design depending on their location in an aluminum smelting plant. 

The situation improved drastically when it was understood that the instability was caused 
by the interaction between the horizontal currents appearing in the aluminum layer in the 
result of the interface deformation and the vertical component of the magnetic held created 
by the external current supply lines js-ll.]. Upon development of effective modeling tools 


(see, e.g., jl2 , 13]) and applying them to designing new and retrofitting existing supply lines, 
the problem was largely solved. 

We should stress that the analogy between an aluminum reduction cell and a liquid 
metal battery is far from complete. Not only a battery has three layers instead of two, 
its aspect ratio, for the laboratory prototypes developed so far, is not small (in fact, the 
optimal geometry of a liquid metal battery is yet to be determined). Furthermore, the typical 
current density is much higher (about 1 A/cm 2 ) in a battery. The analogy is, therefore, just 
a starting point of our analysis. Nevertheless, for the absence of a better name, we will use 
the term ‘metal pad instability’ for the instability mechanisms considered in this paper. 
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The analysis follows the approach, in which one instability mechanism is analyzed sep 


arately fromMhe others (similar approaches were recently used for the Tayler 
convection 




and 


7] instabilities). We also apply a drastic simplification replacing the liquid 
metal layers A and B by slabs of solid metals suspended above and below a liquid layer of 
a poorly conducting electrolyte. The large-scale sloshing motions of the metal layers (for 
example, the gravitational waves) are represented by the motions of the slabs, which we 
model as two-dimensional oscillations of mechanical pendula modified and coupled to each 
other by the electromagnetic forces. The approach is similar to that successfully applied 
to the metal pad instability in the aluminum reduction cells in jlJ ]. We go further than 
simply modifying the results of [l^| to the case of a three-layer system. A broader range 
of possible interactions between the currents caused by the interface deformation and the 
imposed magnetic fields is considered. 


II. MODEL 

A. Simplifying assumptions 

The system analyzed in this work is shown schematically in Fig. [21 The metal layers of 
a battery are represented by solid metal slabs A and B rigidly attached to weightless rigid 
struts pivoted at the top. The free oscillations of the slabs around the two horizontal axes 
passing through the pivot imitate the sloshing motion of the liquid layers. The slabs are 
separated from each other by a layer of a poorly conducting electrolyte. Strong electric 
current of density J flows through the system. There is also the magnetic field of induction 
B , which combines the field induced by J and the field induced by currents in external 
supply lines and neighboring batteries. 

In the rest of the paper, symbols without superscripts are understood as related to both 
pendula. To indicate, where necessary, the properties and variables related to a specific 
pendulum, superscripts A or B are used. 

In the unperturbed state, the slabs’ surfaces are horizontal, the thickness of the electrolyte 
layer is constant h 0l and the electric current is uniform and purely vertical: 

J = J 0 = —Joe z , J 0 = const > 0 (1) 

(for consistency, we will always consider a battery in the process of being charged, but the 
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FIG. 2: Model mechanical system: two independently suspended metal slabs are separated by 
a layer of liquid poorly conducting electrolyte. In the unperturbed state, the electrolyte layer 
has thickness ho, and the constant vertical current Jo hows through the system. Motion of the 
slabs perturbs the thickness by rj(x,y,t) and causes perturbations of electric current j. Pictures 
illustrate, schematically and in an exaggerating way, the results of the tilting of the slabs around 
the y-axis (a) and x-axis ( b ). 

derivations and results are equally valid for a discharging battery with Jo < 0). When 
any of the slabs is tilted, the local thickness and, thus, local resistance of the electrolyte 
changes, and the current perturbations j(x, t ) appear in the slabs and the electrolyte. Their 
interaction with the magnetic field creates Lorentz forces that modify the oscillations of 
the slabs, making them coupled with each other’s and different from the purely gravita¬ 
tional oscillations. The hypothesis we explore here is that this effect may be a source of an 
instability. 

The replacement of liquid metal layers by solid slabs is the replacement of a system 
with infinitely many degrees of freedom by a system with just four such degrees. This can 
be considered as a low-mode approach, in which the key physical mechanism, namely the 
coupling between the deformation of interfaces and the electromagnetic forces, is retained, 
and the large-scale dynamics of the system is analyzed. 

We assume that the electric conductivities of the electrolyte and metal slabs satisfy 

o E < o A ~ o B - ( 2 ) 
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The electrolyte layer is assumed shallow, with its thickness much smaller than the typical 
horizontal size of the slabs 

h 0 <C L. (3) 

The perturbation of the local thickness of the electrolyte layer rj(x,y,t) satisfies 

V < ho. (4) 

The first-order approximation in terms of the perturbation amplitude and of the ratio 
Hq/L is used. 

Perturbations of the magnetic field induced by the current perturbations j are assumed 
much weaker than the base magnetic field B and neglected in the analysis. The diffusion 
effects, such as the viscous friction and Joule heat, as well as the pressure forces arising in 
the electrolyte are also neglected. Finally, we assume that for each pendulum the distance 
between the pivot and the center of mass is equal to /to. 


B. Governing equations 


For each slab, we will use the local Cartesian coordinate system rigidly attached to it 
and having the origin at the center of mass. The z-axis is directed upwards along the strut, 
while the horizontal axes are along the main axes of inertia. The motion is described by the 
angular momentum equations for rotations around the horizontal axes passing through the 
pivot and parallel to x and y at the moment when the slab is in the bottommost position: 


L yy 


d 2 e x 

dt ' 2 
d 2 9y 
dt 2 


7~g,x T 7 ~L,xi 
T 9 ,y + T L,y 


( 5 ) 

( 6 ) 


Here, 9 X and 9 y are the angles of rotation (see Fig. [2]), I xx and I yy are the moments of inertia, 
and the right-hand sides are the sums of the net torques of the gravity (r g>x and r Si3/ ) and 
Lorentz {tl, x and tl, v ) forces with respect to the pivot. 

The moments of inertia are 


I xx = M 


dyy ~ M 


K + 

h 2 n + 


L 2 y + H 2 
12 

Ll + H 2 
12 


L 2 + H 2 

M— -, 

12 ’ 

Ll + H 2 


M 


12 


( 7 ) 

( 8 ) 
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for a rectangular slab and 


r r r „ r / l2 3 R 2 + H 2 \ 3R 2 + H 2 

I XX — lyy — Irr ~ M yh Q H-—- J ~ M -—- (9) 

for a cylindrical one. In these expressions, M is the total mass, L x , L y , and R are the 
horizontal dimensions or radius, and H is the height of the slab. 

The torque of the gravity force is 


Tg, x = ~gMh 0 e x , 

r g ,y = -gMh 0 9y. 

The pure gravitational oscillations have the squared frequencies 

gMh 0 12 gh 0 


l ) 


(^11 ) 


hi 

gMh 0 


L yy 


Ly + H 2 ’ 
I2gh 0 
L 2 X + H 2 


for a rectangular slab and 


for a cylindrical one. 


\^x) — \Uy) — 


gMho 12 gho 


3 R 2 + H 2 


( 10 ) 

( 11 ) 

( 12 ) 

(13) 

(14) 


We will use u A to make the equations non-dimensional. Denoting the non-dimensional 


time as t' = tu A , we obtain: 


•2 qA 


d 2 e 


X 


dt' 2 


+ 


A\ 2 
x_ 1 qA 
A 
y 


et = 


T, 


x,L 


(uj a ) 2 I A ’ 

\ y ) ± xx 


Ml , „A 


T, 


y,L 




2 OB 


d 9 


dt' 2 

d?9y 

~M 2 


+ 


' ,B\ 2 

~“T I = 


T, 


B 


x,L 


id, 


y 


(tu A ) 2 I B ’ 
\ y ) ± xx 


c u: 


+ 1^1 e? = 


T, 


y,L 


UJi 


«) s 


(15) 

(16) 

(17) 

(18) 


In order to compute the torque of the Lorentz force, we need to find the perturbations 
j of the electric current caused by the tilting of the slabs, specify the magnetic held B, 
compute the force density f = j x B, and integrate the components of its torque 


TL, x Vlz^rhofy-) 

(19) 

7~L,y %fz h'ofx* 

(20) 
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The first step of this procedure is discussed here. The rest is completed for the specific 
cases of our analysis in section IIIII 

The local thickness of the electrolyte is, in the asymptotic limit of low-amplitude pertur¬ 
bations (see Fig. [2|), 

h = h 0 + 7 i(x, y, t) = h 0 + ( 9y - 0 A ) x - (0f - d A ) y. (21) 

We can always choose the axes so that, at a given moment of time, the thickness is given by 

h = h 0 + rj(x, t) = h 0 + (0f - 0 A ) x. (22) 

If the Lorentz force has zero torque component r EtX in such coordinates, the subsequent 
oscillations of the slabs occur around the y- axis, i.e. so that (122]) remains valid. 

The variation of the electrolyte thickness causes variation of its local electrical resistivity 
and, thus, perturbations of electric currents. To evaluate them, we use the assumption ([2]) 
and require that the surfaces of the slabs facing the electrolyte remain equipotential: 

<f> = 0 at B, <f> = <f> A = const at A. (23) 

Furthermore, as a first-order approximation, we assume that the value of does not change 
when the interface is tilted: 

<f> A » (24) 

and that the perturbed currents in the electrolyte layers remain vertical: 

J E = —JqZz + jf(x, t)e z . (25) 


Considering that 


Jn = 


%aE and J„ 


3? = 


& A cr E 


and using 


ho ho + rj 

, we find the distribution of current perturbations in the electrolyte: 


■ E §oO E $n <7e . . $o '°E 


Jz = 


h 0 + r) 


h 2 

n o 


V = Jc 


V 


( 26 ) 


We now derive the expressions for the current perturbations within the solid slabs. The 
derivation is first conducted for the bottom slab B. We employ the fact that, since the 
electric conductivity is high, the current perturbations can be assumed to form completely 
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closed loops within the slabs. This implies the boundary conditions: 


j B I = jE 

Jz \ z= J{B /2 Jz 1 

1 B \ = 0 

Jz \ z = _ H B /2 

it ■ n L = °> 


(27) 

(28) 
(29) 


where —H B J 2 < z < H B / 2 is the vertical coordinate within the slab, j B = (j B ,j B ), and 
is the normal to the slab’s boundary <9f2 in the x-y- plane. 

The vertical component is approximated as 


n 


■B 


-n E = T^L 

2 3z J °2 h 0 ' 


(30) 


The derivation of the horizontal currents j B uses the condition of zero free charges 


V-j B = 0 


and the vertical integration 


r H 3 / 2 

~jl = / Hdz. 

J-H 3 / 2 


Integrating (T3T1) and applying (127]) - (1251) . we find 


= _,b 

JZ * 


+ 


dx dy 

In the simpler case when the oscillations occur in one plane 
bation currents have only the x-component, and (1331) integrates to 

Jo 


3x = 


2hr 


(d B - dy) x 2 + const. 


(31) 

(32) 

(33) 

the horizontal pertur- 

(34) 


In the interesting for us case of cylindrical slabs (see section IIII BD . the boundary condition 
(1291) leads to 

if = ^ W - K) (R 2 - r 2 ) . (35) 

where r = ( x 2 + y 2 ) 1 ^ ■ 

In the general case of two-dimensional oscillations (121]) . the derivation is slightly more 
complex. We will need the currents in rectangular slabs in section IIII A1 Integrating (l33|) in 
the ^-direction and using (129]) . we find 


dll 

dx 


r L v / 2 
l-Ly/2 


qB _ qA 

j z dy = — y — —- L y J 0 x , 


h. 


(36) 
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where 


rLy/2 

!x(x) = Jxdy 

J-Ly/2 

is the y- 2 -integrated (E-component of the current. Integrating (135j) along x and applying 
again we find 


I B = 

X 


LyJo ( a B 
2 ho 


&-»}) 


X 


( 37 ) 


In the same manner, we obtain the distribution of the x-z-integrated y-conrponent 

2 


I B = 
y 


L x Jq / q b 


2 hn 


(sf - <£) 


T ] - y 2 


(38) 


Following a similar procedure or simply applying the charge conservation condition, we 
find the currents in the upper slab A: 


a a _ _ -n tA _tB jA 

J± J ±1 1 x 1 x 1 1 y 


-I B j A = j 

y ? Jz J. 


B 


(39) 


This completes the preparatory derivations. 


III. SOLUTION 


To complete the governing equations and start solving the problem we need to specify 
the magnetic field B(x). In a real battery, B is a complex three-dimensional field, which 
includes the component induced by the base current J 0 and the components induced by 
the currents in the electric supply lines and, if present, neighboring batteries. Since we 
solve a linear problem, the analysis can be simplified and given clearer physical meaning 
by conducting it separately for selected components of B. We start, in section IIII Al with 
the interaction of a purely vertical magnetic field and horizontal currents, i.e., with an 
analog of the solution 14] for the mechanical model of an aluminum reduction cell. Section 
IIIIBI presents the more interesting results dealing with the interaction between the current 
perturbations and the azimuthal magnetic field induced by J 0 . 


A. Case 1: Vertical magnetic field 

The driving mechanism of the metal pad instability in the aluminum reduction cells is 
the interaction between the horizontal current perturbations and the vertical component of 
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the externally (by the neighboring cells and supply lines) generated magnetic field BQ- 
The interaction creates electromagnetic coupling between the gravitational waves at the 
aluminum-electrolyte interface and causes the instability. 

In order to explore the possibility of an analogous instability in a liquid metal battery, 
we consider a system with rectangular metal slabs and assume the imposed magnetic field 
of the form 

B = B 0 e z , B 0 = const. (40) 

Taking the cross-product with the integrated currents f[3Tl) . (j38lh we find, for the slab B, 
distributions of the correspondingly integrated Lorentz force components along the x- and 
y- axes: 

F*{y) = Iy(y)B 0 , F B (x) = -I B (x)B q (41) 

and of the torque 


Tl,x( x ) = h oFy = - B ° J ° Ly (Q b - 0 A ) 


tlM = ~h 0 F x = 




— x 


- y 


(42) 


(43) 


Integration along the respective coordinates gives the final expressions for the net torque: 


r L = 


r B = 
T L,y 


rL x /2 

' —Lx/2 

r-Ly/2 


T L)X {x)dx = 


BqJqL v L^ 


J-Ly/2 

For the slab A, we have (see (J39]) ): 


r L , y (y)dy = 


12 

QjQlJx-^y 


w-<a 


BnJoLrL 3 


12 


(<£ - £) • 


(44) 

(45) 


A = —t B t A = —t B 
'L,x ‘L,xi 'L,y ‘L,y’ 


(46) 


Substituting QS1) - Q46]) into the non-dimensional governing equations (fT5]) - (lT8l) . we obtain 

d 2 9 A ' 4x2 


dt’ 2 


d 2 e B 

dt 12 

d 2 e B 
dt 12 


+ (j^) 0 . A = « - Of ). 


d 2 e A 
~di i2r 

,B\ 2 


+ 9* = - e A Gi (0 B - 9*) , 


+ (^) ~ ~ eBG * ( e v ~ e v ) ’ 

nB C B r* B (ft B fi A \ 

Vy C '- J y (y x v x J , 


(47) 

(48) 

(49) 

(50) 
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where 


s~iA 

— 




G? = 


G? = 


lI + (h a y 


L 2 y + (H A ) 

S 

LV 

li + (h a y 


2 ’ 


l; + (h b ) 


2 ’ 


r2 t 2 


+ (H A Y 


L l L\ + (H B )‘ 


are the non-dimensional geometry factors, and 


e A = 


BMl 


12p A gh 0 H A ’ 


e s = 


BMl 


12 p B gh 0 H B 


(51) 

(52) 

(53) 

(54) 

(55) 


are the non-dimensional control parameters that determine the strength of the electromag¬ 
netic effect. 

The solution of f|47p ~ f|50p is a linear combination of the eigenmodes 


6 ~ exp ( 77 1 ') 


(56) 


with eigenvalues 7 = u±ia. The real part u j is an electromagnetically modified gravitational 
frequency. The imaginary part < 7 , when non-zero, is the growth rate of this eigenmode. Pres¬ 
ence of at least one pair of complex-conjugate eigenvalues in the spectrum implies instability. 

It is convenient to further simplify the geometry of the system, so that we can obtain 
an analytical solution of the problem. One possibility is to consider a battery, in which 


H A = H B and, thus, u A = lo b = c u x , u A = u B = u y , G A = G A = G x , G A = G B = G. 


yA 

T x 


-M 

T x 


Introducing the new variables 


& = - C P, = Oy - 0„ 


(57) 


and subtracting (1471) from (l49li and (|48l) from (1501) we obtain the reduced system 

d 2 (3 x / x2 


dt' 2 


+ Px = ~G X (t A + e B ) Py, 


d 2 P y 

dt' 2 


— Py — Gy (e* 4 + 0 x . 


(58) 

(59) 

(60) 
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Substituting (1561) and solving the quadratic equation for q 2 we find positive determinant 
and roots q 2 > 0 if 2 ( e A + e B ) ( G x G y ) 1 ^ 2 < |l — (u x ) 2 / (u y ) 2 \. The system is stable in this 
case. On the contrary, if 


2(e A + e B ) (G x G y ) 1/2 


> 


-1 ^X 


(61) 


the solution necessarily has a pair of complex-conjugate eigenvalues q = u ±ia, a > 0, and, 
therefore, the system is unstable. 

In order to relate our results to those found for the aluminum reduction cell, we take the 
same asymptotic Hmit as m Q. We assume that uot only h 0 , but also the thicknesses of 
the metal slabs II ', II 11 are much smaller than the horizontal dimensions L x and L y . This 
leads to 

12 h 0 g 


, s 2 12 h 0 g 


l Uy) = 


LI 


T 2 L 2 

_ _ y_ 

~ i ,2 — 1,2 


(62) 


+ 


BqJq 


> -w*|. 


(63) 


by the presence of two terms in the 


and allows us to rearrange (j6TT) as 

2 ( BoJ* 

\p A H A ' p B H B i 

which only differs from the instability criterion in |l^| 
left-hand side. 

In the general case, the system (1TH) - (|501) is not reducible to two equations and does not 
have a simple analytical solution. Considering, however, that the ratios of the gravitational 
frequencies and the geometric factors are all of the order one, we can write the instability 
criterion, approximately, as 


C A e A + > 


B 


1 - 


OJZ 


ur, 


(64) 


where C A ~ 1 and C B ~ 1 are the constants accounting for the effect of geometry. 

The physical interpretation of the criterion is similar to the interpretation of the metal 
pad instability in aluminum reduction cells [blj]. The instability occurs in a battery when 
the product JqBq exceeds a limit determined by the geometry of the battery. The limit 
decreases linearly with decreasing thickness of each layer: ho, H A , or H B . It also depends 
on the horizontal shape of the battery. A battery of square cross-section L x = L y and with 
H A = H B has u x = u y and, so, is always unstable (the same conclusion can be easily shown 
as valid for a cylindrical battery). In general, we expect that, among the batteries with 
given H a and H B , the square and cylindrical ones would be most prone to the instability. 
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The critical values of e A and e B would increase with increasing difference between L x and 

Ly. 

An essential difference between the instabilities in the aluminum reduction cell and the 
battery is manifested by the combination of e A and e B appearing in (1(3-0) in place of just one 
such parameter. The presence of the second metal layer makes the system more unstable. 

B. Case 2: Azimuthal magnetic field 

Liquid metal batteries are different from aluminum reduction cells in many respects other 
than the presence of the top metal layer. In particular, the density of the base electric current 
Jo is about an order of magnitude higher, and the ratio between the vertical and horizontal 
dimensions of the metal layers is not necessarily small. This alters the electromagnetic 
interactions and may activate new mechanisms of instability. In this section, we demonstrate 
such a mechanism. 

The instability is caused by the interaction between the current perturbations and the 
azimuthal magnetic field induced by Jo- For simplicity, we consider a cylindrical battery, in 
which the magnetic field is 

B 0 = Boe# = = - sin <pB 0 e x + cos (t>B 0 e y , (65) 

where /io is the magnetic permeability of free space, and (j) is the polar angle in the x — y- 
plane. As will be seen from the following discussion, a similar instability should appear in 
a battery of an arbitrary cross-section. 

The coordinate system is oriented so that the deformation of the electrolyte thickness at 
some moment of time is along the x-axis, i.e., described by (|22j) . The current perturbations 
in each slab have only two components j x and j z (see (l30j) and ([35]) ). The instantaneous 
distributions of the Lorentz forces integrated over the thickness of each slab are: 

f A = j A x B 0 = [H A x 2 e x + H A xye y + x (R 2 - r 2 ) e z \ ( 9 y - 9 A ) , (66) 

f B = j B x B 0 = [H B x 2 e x + H B xye y - x (R 2 - r 2 ) e z ] (i 9 B - 9 A ) . (67) 
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The components of the torque are: 

= '~4K ^ _ T ' 1 ) + HAh ° xyl i W ~ °v) • < 68 ) 

ft = ^ h* 2 ( R 2 - r 2 ) - ff-%* 2 ] (»* - e$) , (69) 

if* = [~ x y (-R 2 - r 2 ) + H B h„xy] (0 B - 0*) , (70) 

fly = [- 2 (« 2 - r 2 ) - ff%* 2 ] (Of - «?) ■ (71) 

Integrating them over the slab, we hnd 



The torque of the gravity force has the components (see (TlOT) — (TlTh ): 

r g,x = T tv = °> T iv = - gMh 0 6 ,r£ y = -gMh 0 9 (76) 


We see that neither slab experiences torque around the x-axis. The oscillations will 
remain in the x — z-plane and the governing equations are reduced to those for just two 
degrees of freedom: 0 y and 6 y . The non-dimensional equations (fT5]) -(fT8 ]) can be rewritten 
as 

^ + K = H" - K ") K ~ <) • < 77 ) 

r l'2nB 

_i + G «B (78) 


where we have introduced the non-dimensional geometry parameter 



and the non-dimensional control parameters 

a _ ^qJqR 4 

6 ~ 48 p A ghlH A ' 

A _ Mo^o-R 2 

" ~~ l6p A gh 0 ’ 


B _ /i() J‘o -R 4 

6 = 48 p B gh 2 0 H B ' 

B _ A t O'7o-R 2 

K ~ 16p s A 


(79) 


(80) 

(81) 
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that evaluate the strength of the electromagnetic torque caused by the horizontal (j80|l and 
vertical (l8Tf perturbations of currents. Since ho is much smaller than R and H , it is safe to 
assume that e A > n A and e B > k b . 

The rest of the analysis is straightforward. We use the ansatz 9 A , 9 B ~ exp t'yt' and solve 
the quadratic equation for y 2 . The determinant of the equation is always positive: 

D = (s A + s B ) 2 - 4s A s B + 4Ga A a B > 0, (82) 

where a A = e A + n A , a B = e B — k b , s a = 1 — cv 4 , and s B = G{1 — a B ). Of the two roots 

7 h = \ O + ^ ± -0 1/2 ), (83) 

the larger is always positive and corresponds to the oscillations of the slabs with the frequen¬ 
cies modified by the Lorentz forces. The smaller root becomes negative if Ga A a B > s A s B , 
which can be rewritten as 

e A + n A + e B -k b > 1. (84) 

Since such a root corresponds to the presence of an eigenmode 9 A ,9 B ~ e at ' with a = 
(—y 2 ) 1 ^ 2 > 0, the condition (154)1 is a criterion for instability. According to it, the effect 
of the horizontal current perturbations is always destabilizing. The effect of the vertical 
current perturbations is weaker and mixed: destabilizing for the slab A and stabilizing for 
the slab B. 

IV. DISCUSSION 

The predictions made in this paper are, strictly speaking, valid only in the framework of 
our mechanical model. At the same time, they are expected to be qualitatively applicable 
to operation of real batteries. This includes the principal physical mechanisms of the insta¬ 
bilities, the form of the control parameters (T55|) (j80lh and (15T]h and, possibly, the order of 
magnitude of the critical values of these parameters. The results of this paper can serve as 
a stating point for future work based on more realistic models. 

By analogy with the aluminum reduction cells, we predict that the instability of the type 
described in section IIII Al will occur in real batteries if sufficiently strong vertical magnetic 
fields are allowed. The batteries of square or round horizontal cross-section will be particu¬ 
larly unstable. A more specific prediction requires detailed analysis of specific geometries. 
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We can list three factors that make the instability in the batteries more likely than in the 
reduction cells. One is the about an order of magnitude higher density of the base electric 
current Jo- Another is the presence of the top metal layer, which plays a destabilizing role, 
and whose parameter e A is particularly large because of the low density p A . Finally, thin (a 
few mm) electrolyte layers are used in the current battery prototypes. This can be compared 
with the thickness about 4 cm in modern reduction cells. 

At the same time, the horizontal dimensions of a battery are unlikely to be as larger 
(several meters) as in typical reduction cells. Finally, the instability can be postponed or 
even completely avoided by optimizing the current supply lines so that there is no significant 
vertical magnetic held. 

We have to be more careful while making predictions for the instability described in 
section IIII Bl On one hand, it appears to be more dangerous than the instability of the first 
type, since the azimuthal magnetic held (165]) cannot be ‘optimized away’. On the other 
hand, existence of such an instability in a real battery yet needs to be confirmed. 

At this point, we assume that the instability of the second type appears in real batteries 
and make preliminary estimates of the size, at which this would happen. Since ho <C 
R, H a , H b , and p A is about an order of magnitude smaller than p B , e A is much larger than 
the other parameters in (l84]h For simplicity, we approximate the instability criterion as 
e A > 1. Using J 0 = 10 4 A m -2 and p A = 500 kg m~ 3 (approximate value for liquid lithium 
at 720 K), we find the critical values of radius, above which the battery is unstable, shown 
in table Q] They are comparable with the typical critical radii predicted for the Tayler m 
or convection [7] instabilities. 


ho 

[mm] 

1 

1 

5 

5 

H 

a /R 

1 

1/4 

1 

1/4 

R cr [m] 

0.12 

0.078 

0.36 

0.227 


TABLE I: Critical radius, above which the battery is unstable to the instability caused by the 
azimuthal magnetic field. See text for explanation. 

Our final remark is that the metal pad instabilities may play a substantial role in the 
operation of scaled-up liquid metal batteries and have to be included into the future analysis. 
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